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LOW  ALTITUDE  PLASMA  LINE  ANISOTROPY 


INTRODUCTION 

Plasma  line  observations  obtained  from  incoherent  radar  backscatter  have  been  used  as  a 
ground-based  method  for  deriving  information  about  the  size  and  anisotropy  of  the  ionospheric 
photoelectron  fluxes.  In  the  past,  interpretation  of  plasma  line  data  has  been  confined  to  alti- 
tudes above  the  F2  peak  where  plasma  lines  enhanced  by  photoelectrons  moving  upward  have 
larger  amplitudes  than  those  enhanced  by  downward  moving  photoelectrons.  The  net  upward 
electron  flux  in  the  topside  ionosphere  arises  because  the  photoelectron  mean  free  path  is  large 
and  within  a mean  free  path,  more  electrons  are  produced  just  below  a given  point  of  observa- 
tion than  just  above  it. 

At  lower  altitudes  the  photoelectron  mean  free  path  is  small  compared  to  the  gradient 
scale  length  of  the  photoelectron  production  rate,  yet  plasma  line  anisotropy  is  still  observed. 
Measurements  below  the  peak  were  first  shown  by  Yngvesson  and  Perkins  (1968)  and  recently 
by  Legeune  and  Kofman  (1977).  These- observations  have  consistently  shown  anisotropy  in 
the  plasma  line  amplitude  ranging  from  20%  to  a factor  of  two  for  upgoing  to  downgoing  parti- 
cles. Measurements  made  at  Arecibo  have  shown  a 20-50%  anisotropy. 

In  this  paper  we  describe  calculations  of  the  plasma  line  intensity  which  use  a detailed 
transport  equation  solution  for  the  photoelectron  distribution  function.  We  postulate  from  the 
work  described  below  that  the  observed  low  altitude  anisotropy  reflects  the  forward  peaked 
property  of  the  elastic  electron  neutral  scattering  cross  sections.  We  further  argue  that  given 
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anisotropic  cross  sections,  an  anisotropy  in  the  low  altitude  electron  fluxes  arises  when  the  lo- 
cal mean  free  path  is  of  the  order  of  the  local  scale  height. 


LOW  ALTITUDE  FLUX  ANISOTROPY 

The  photoelectron  transport  calculation  used  in  our  plasma  line  studies  has  been 
described  in  Oran  and  Strickland  (1976).  Further  technical  description  of  the  general  method 
can  be  found  in  Strickland  et  at.  (1976).  Figures  1 and  2 show  results  of  a typical  calculation  of 
the  photoelectron  flux  as  a function  of  pitch  angle  at  157  and  3SS  km  for  selected  energies. 
Figure  3 shows  the  way  in  which  the  calculated  flux  anisotropy  of  4 ev  electrons  varies  with  al- 
titude. 

The  low  altitude,  the  low  energy  flux  anisotropy  shown  in  Figures  1 and  3 arises  as  a 
consequence  of  the  fact  that  the  mean  free  path  of  an  energetic  photoelectron  (E  > 10  ev)  can 
be  of  the  same  order  of  magnitude  as  the  local  neutral  density  scale  height.  This  fact,  coupled 
with  an  anisotropy  in  the  elastic  electron-neutral  collision  cross  sections,  results  in  a net  down- 
ward flow  of  particles.  The  cross  sections  for  atomic  oxygen,  which  are  shown  in  Figure  4 
(Blaha  and  Davis,  1975),  show  a typical  angular  dependence;  they  are  forward  peaked  and 
reach  a minimum  at  some  angle  greater  than  ninety  degrees. 

Some  insight  into  the  nature  of  this  low  altitude  flux  anisotropy  can  be  obtained  by  look- 
ing at  a simplified  photoelectron  calculation  based  on  the  limiting  case  in  which  electron- 
neutral  collisions  arc  assumed  to  be  entirely  forward  scattered.  We  lake  the  photoelectron 
equation  of  motion  to  be 

V «'^„(2)v 

(1) 

where  is  the  electron-neutral  collision  frequency  and  z is  the  spatial  coordinate  in  the  verti- 
cal direction.  Then  we  consider  a neutral  gas  in  one  dimension  with  a density  gradient  such 
that 


2 


•Jl 

i 

i 


I 


Vf„(Z)  •vge 

where  L is  the  neutral  density  scale  height.  We  then  wish  to  solve  an  equation  of  the  form 

BF.  BF. 

dz  av 

which  represents  a steady  state  one  dimensional  Boltzmann  equation  for  the  distribution  func- 
tion, Fp,  with  a source  of  electrons  at  velocity  ±Vg.  We  assume  that  the  electrons  move 
through  a medium  in  which  their  mean  free  path  is  very  small  compared  to  the  gradient  scale 
length  of  their  production  rate.  This  assumption  is  valid  in  the  atmosphere  up  to  about  200 
km  and  allows  us  for  the  purposes  of  this  discussion  to  set  equal  to  a constant 


The  solution  to  Eq.  (3)  is 


(4) 


where  t>g  is  the  collision  frequency  at  z —O.  The  superscripts  on  /^,  -I-  and  — , denote  upgoing 
and  downgoing  particles,  respectively.  With  no  loss  of  generality  we  can  define  the  altitude  of 
interest  to  be  z — 0,  and  Eq.  (4)  becomes 


which  is  valid  when  (v„L  ~v^)  >0.  This  condition  guarantees  that  electrons  created  at 
2—00  cannot  reach  z — 0.  Violation  of  this  constraint  leads  to  a singularity  in  f~,  since  /3 
does  not  vanish  at  «>  in  this  simple  model.  We  define  an  anisotropy  parameter 


Jfi  n **0  + - lvl 

”//(|v|)  “ % - -|v|)/L- 

Note  that  there  is  no  anisotropy  when  |v|  ~ Vg  since  A (Vg)  » I.  The  anisotropy  builds  up 
as  the  electron  energy  degrades. 
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In  order  to  concentrate  specifically  on  the  low  energy  particles,  we  define  the  local  mean 
free  path  as  = \glv„,  so  that  when  Z,'^'7L  < 1, 


I . . . 1 + L'^'PIL 

V <<  V-)  = 

1 -L'^P/L 


(7) 


which  shows  that  a significant  anisotropy  exists  at  low  energies  when  the  local  photoelectron 
mean  free  path  is  comparable  to  the  scale  height.  For  example,  from  Figure  5 we  see  that  this 
condition  is  satisfied  for  10  ev  electrons  at  ISO  km. 


We  emphasize  that  Equation  7 is  the  result  of  a qualitative  calculation  given  to  help 
understand  the  detailed  simulation  results  shown  in  Figures  1,  2,  and  3.  Aside  from  the  un- 
realistic source  term.  Equation  3 includes  only  forward  scattering  with  no  diffusion  effects. 
Thus  it  overemphasizes  the  anisotropy.  The  detailed  simulation  is  a solution  of  the  steady  state 
Boltzmann  equation 


(8) 


where  0 is  a local  source  term,  v represents  the  loss  of  energy  by  the  photoelectrons  to 
the  ambient  electrons  through  collisions  and  plasma  wave  excitations,  and  dF  / d/j^  involves 
the  details  of  elastic  and  inelastic  collisions  among  electrons  and  neutral  species.  Solution  of 
Equation  (8)  involves  following  the  detailed  degradation  of  the  primary  electron  spectrum  pro- 
duced by  ionization  due  to  solar  radiation  and  yields  the  distribution  function,  Fiz.vji)  where 
M is  the  pitch  angle  (Oran  and  Strickland,  1976). 


PLASMA  LINES  FROM  INCOHERENT  RADAR  SCATTER 


When  the  radar  wavelength,  X,  is  very  much  larger  than  the  Debye  length,  the  radar  sig- 
nal is  scattered  incoherently  by  electron  density  fluctuations.  Because  the  ionosphere  is  a 
non-equilibrium  plasma,  the  intensity  of  the  electron  density  fluctuations  is  enhanced  beyond 
that  expected  from  purely  thermal  fluctuations.  The  fast  photoelectrons  feed  energy  into  and 
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excite  the  plasma  waves  which  in  term  are  Landau  damped.  The  total  amplitude  of  the  plasma 
line  signal  represents  an  equilibrium  between  a number  of  processes  which  excite  and  damp 
the  plasma  waves  scattering  the  radar  beam.  Collisions  between  thermal  electrons  and  ions 
can  excite  plasma  waves  and  the  interaction  with  energetic  photoelectrons  can  enhance  the  sig- 
nal amplitude  above  the  thermal  equilibrium  value.  The  waves  are  damped  by  electron  colli- 
sions with  ions  and  neutrals  as  well  as  by  Landau  damping. 

The  plasma  line  signal  amplitude  can  be  expressed  in  terms  of  a temperature,  Tp,  which 
reduces  to  T,  when  the  plasma  is  in  thermal  equilibrium.  The  upshifted  (downshifted)  signal 
is  a result  of  scattering  off  of  plasma  waves  excited  by  electrons  moving  primarily  down  (up) 
the  magnetic  field  line.  Thus  the  returned  signals,  proportional  to  and  Tp~ , may  differ  by 
an  amount  related  to  the  degree  of  anisotropy  in  the  photoelectron  flux  spectrum.  The  balance 
of  processes  exciting  and  damping  the  plasma  waves  leads  to  an  expression  for  Tp, 

_ X"  +fp  +f„t 

^ + f„b  + Lp  (9) 

(Perkins  and  Salpeter,  1965a  and  b).  In  this  equation  and  fp  are  the  Maxwellian  and  pho- 
toelectron contributions  to  the  plasma  line  source  term  and  represents  the  effects  of  elec- 
tron collisions  with  ions  and  neutrals.  The  scattering  occurs  off  of  plasma  oscillations  of 
wavelength  X/2  and  phase  velocity  v^  and  “y  The  returned  plasma 

line  signal  is  Doppler  shifted  about  the  radar  frequency  by  an  amount  ±vr,  where 

, , l2kT, 

“ •'p  + ~2 + sine  = 

^ '"e  (10) 

where  pp  is  the  plasma  frequency,  nip  is  the  electron  mass,  t>„  is  the  electron  cyclotron  fre- 
quency, and  9 is  the  angle  between  the  direction  of  radar  beam  propagation  and  the  magnetic 
field. 
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In  the  presence  of  the  earth's  magnetic  held. 


^ ^ I 1/2  oo 

^1  2:  ^ I„(bsin^0)  e -(>--«) V2cos^« 


e } fl  *■  —00 


(Selpetcr,  1961),  where  Ng  is  the  ambient  electron  density,  1„  is  the  Bessel  function  of  ima- 


ginary argument. 


y ~>'R/‘'ce- 

_ 4 

" khl' 


n “ —00  '^min 

(Fremouw  et  a!.,  1969),  where  J„  is  the  Bessel  function  of  real  argument, 

2 V,  sin  f) 

X*/  _ ^ 


2 cos » ■ 

The  Vi  and  u are  the  components  of  electron  velocity  perpendicular  and  parallel  to  the  magnet- 
ic field,  and  F{\^,u)  is  the  three-dimensional  velocity  distribution  function  for  electrons  where 
w is  evaluated  at  «„.  The  lower  unit  of  the  integral,  v^jp,  corresponds  to  that  velocity  at  which 
the  photoelectron  distribution  function  begins  to  deviate  substantially  from  Maxwellian.  The 
photoelectron  contribution  to  the  damping,  L^,  can  be  expressed  as 

2-nkT,  “ /•»>  . II  ( ^ a _ a 1 \ 


y cos 


**  n-~oo ^ mgU  du  mg\^  dv^^J  ^ u-u. 


Substitution  of  a Maxwellian  distribution  function  into  and  both  yield /,„j.  The  x”  in 
Equation  (9)  is  related  to  the  electron-ion  and  electron-neutral  collision  frequencies  by 
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(19) 


where  Dis  the  Debye  length. 

The  plasma  line  is  in  thermal  equilibirum,  Tp  — when  either  thermal  electron  Lan- 
dau damping  dominates  or  collisional  damping  dominates.  Thus  there  is  only  a limited  range  of 
phase  energies  for  which  photoelectrons  are  the  dominant  contribution  to  the  plasma  line. 

Figure  6 summarizes  a set  of  plasma  line  data  taken  at  Arecibo  Observatory  in  February 
27,  1977.  The  Arecibo  Observatory  is  part  of  the  National  Astronomy  and  Ionosphere  Center, 
which  is  operated  by  Cornell  University  under  contract  with  the  National  Science  Foundation. 
The  plasma  line  measurements  were  interleaved  with  measurement  of  electron  and  ion  tem- 
peratures and  electron  densities.  We  note  that  the  anisotropy  in  plasma  line  temperatures 
shown  in  this  data  is  typical  of  the  20-50%  enhancement  of  7’^'*’  over  Tp  found  at  Arecibo. 
Uncertainties  in  the  data  points  are  small  and  arise  only  through  calibration  uncertainties.  During  tlie 
nighttime,  the  asymmetry  in  the  upshifted  and  downshifted  plasma  lines  disappears.  Superimposed  on 
the  figure  are  the  measured  values  of  Vjf. 

RESULTS 

The  approach  we  have  taken  is  to  calculate  the  magnitude  and  pitch  angle  anisotropy  of 
the  photoelectron  distribution  function  from  the  detailed  transport  equation  and  then  derive 
the  plasma  line  temperature  we  would  observe.  In  this  way  we  do  not  expect  the  calculation  to 
be  restricted  to  high  altitudes  or  a particular  range  of  phase  energies.  Figure  (7)  shows  calcula- 
tions of  Tp'*'  and  Tp~  made  with  the  photoelectron  distribution  function,  F(z,\,ti),  and  Equa- 
tions (9)  through  (20).  The  values  of  superimposed  in  the  profiles  were  determined  by  elec- 
tron temperatures  and  densities  measured  at  Arecibo  Observatory. 
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Figure  (8)  graphically  summarizes  the  important  terms  contributing  to  Equation  (7). 
Below  about  200  km,  f~  is  greater  than  indicating  a small  anisotropy  in  electron  flux; 
more  electrons  are  going  down  than  up.  However  in  this  region  T^,  is  determined  primarily  by 
a balance  between  the  source  term./^,  and  the  loss  term,  Lp.  From  Figure  (8)  we  see  that  up 
to  about  200  km, 


f. 


p 


" /- 

) Jp 


indicating  that  the  relative  difference  between  and  Lp  is  greater  than  that  between 
and  f~.  In  Figure  (9)  we  have  rewritten  the  terms  in  the  denominator  of  Equation  (9)  in 
terms  of  damping  rates. 


At  high  altitudes,  is  considerably  larger  than  f~  and  these  remain  the  important 
source  terms  until  about  400  km.  The  increasing  importance  of  the  damping  due  to  thermal 
electrons  tends  to  reduce  the  importance  of  the  difference  between  Lp  and  L~.  Up  to  that  al- 
titude where  thermal  electrons  are  dominant  and  when  is  very  low,  is  greater  than 


An  important  point  we  wish  to  make  is  that  previous  calculations  of  the  source  and  loss 
terms  due  to  photoelectrons  treat  photoelectrons  with  velocities  less  than  as  Maxwellian. 
We  note  that  especially  at  low  altitudes  the  calculated  photoelectron  distribution  function  devi- 
ates strongly  from  Maxwellian  at  energies  of  a few  ev.  At  altitudes  below  200  km,  our  calcula- 
tions indicate  that  50-75%  of  Lp  and  25-50%  of  fp  is  evaluated  in  the  integral  range  v^j„  to  v,^. 

CONCLUSION 

The  investigations  we  have  made  of  the  photoelectron  flux  in  the  ionosphere  have  un- 
covered an  anisotropy  in  the  low  altitude,  low  energy  electron  flux  spectra.  This  arises  through 
the  conversion  of  a spatial  inhomogeneity  in  the  neutral  atmosphere  into  a velocity  distribution 
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anisotropy.  The  key  to  this  is  the  forward  peaked  property  of  the  elastic  electron-neutral  cross 
sections.  We  conclude  that  the  downgoing  electron  flux  is  slightly  larger  than  the  upgoing  flux 
in  the  altitude  range  from  about  130  to  200  km. 

In  the  altitude  range  described  above,  the  plasma  line  amplitude  is  a balance  between  a 
source  term  which  excites  plasma  waves  to  characteristic  temperatures  appropriate  for  pho- 
toelectrons, and  a Landau  damping  term  L^.  Our  calculations  show  that  the  low  energy  pho- 
toelectrons (£<£^)  play  a significant  role  in  the  damping.  Treating  them  as  Maxwellian 
severely  underestimates  their  density  and  neglects  their  anisotropy. 
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1 ~ Calciilatod  photoclactron  flux  m « function 
of  n,  the  cosine  of  the  pitch  angle,  for  selected  ener- 
gies at  167  km. 


Fig.  3 — CalcttlaUd  photoelectron  flux  at  4 ev  a func- 
tion of  pitch  angle  at  aelected  altitudea. 


H (COSINE  OF  PITCH  ANGLE ) 


Pig.  4 — Differential  croaa  aectioni  for  elastic  scattering  of 
electrons  on  the  ground  state  of  atomic  oxygen  (Blaha 
and  Davis,  197S).  llte  labels  represent  enei^es  in  ev. 


